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Abstract. This article describes, in elementary terms, a generic approach to produce dis- 
crete random tilings and similar random structures by using point process theory. The 
standard Voronoi and Delone tilings can be constructed in this way. For this purpose, con- 
vex polytopes are replaced by their vertex sets. Three explicit constructions are given to 
illustrate the concept. 



1. Introduction 

Apart from symmetry and long-range order, also randomness is needed to provide appro- 
priate structure models in physics or chemistry. A well-known example for an application of 
random points respectively point processes is an ideal gas, which, at each instance of time, 
may be described by a Poisson point process (cf. [9] p. 449f). 

There are also some well-known random tilings with applications in crystallography and 
material sciences, such as the Poisson Voronoi and the Poisson Delone tessellations ([20, 17]). 
A more general description of random tilings in the context of quasicrystals can be found 
in [18, 10]. Recent work of Gummelt [7, 8] is also concerned with establishing a connection 
between randomness and quasicrystalline structures. Most approaches start from randomly 
generated points in Euclidean space IR'^ and then distill the information about the tiles, often 
convex polytopes, out of the positions of the points. (See Figure 1 for an illustration of this 
widely used concept.) 

The theory presented in this article is partly based on recent work of Zessin [22], where 
random tiles are also extracted from random point configurations. However, this approach 
replaces the concept of convex polytopes by discrete, finite subsets of IR"^, called clusters. 
Nevertheless, tilings of convex polytopes can be embedded into this theory by identifying the 
vertex set of a convex polytope with the polytope itself (see Figure 2). This point of view 
might describe the underlying structure of the atoms and molecules of certain materials more 
accurately. In this article the desired local and global properties of the tiles, in our case the 
clusters, enter by means of so called cluster properties. In particular, the first two examples 
given, are to illustrate this concept. Also, although the examples given in this article are 
tilings, the presented theory is not restricted to these. 

While the typical point processes like the Poisson point process with constant intensity 
might be 'too random' to describe condensed matter, this paper gives an explicit example 
construction to go over to a random tiling 'close' to a quasicrystalline one. 
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Figure 1. A Voronoi tessellation: The (random) points generate their 'neighbourhoods' 




Figure 2. We identify a convex polytope with its vertex set. Here, {X) 
denotes the convex hull of the set X = {ai, . . . , 05} C |R<^ 

2. Random Points 

Point processes can be constructed in various spaces. Here, we restrict ourselves to random 
points in R'^. See [12] for a more detailed and more general description. To speak about 
randomness and probabilities, we first need to identify the objects which should be realized 
randomly. Then we need a notion of which events can be computed or, more precisely, can 
be given a probability. 

Technically, the space for the point configurations is the set of the locally finite simple 
counting measures in IR*^, denoted by ^'(R'^). All possible events are subsumed under the 
cr-algebra generated by the counting functions. (A very readable introduction to probability 
theory and thereby an explanation for the need of terms like measurability and cr-algebras 
is given in [6].) A probability measure on ^*(R'^) equipped with this cr-algebra is called 
(simple) point process. 

In essence, we can say that the point configurations we want to consider have only finitely 
many points in every bounded subset of IR^, and that the typical events are of the form 

I?? G ^•(R'^) : rihask points in , (1) 

where B is some bounded subset of IR*^ and A; is a non-negative integer. A point process may be 
viewed as a mechanism to randomly generate point configurations obeying a given probability 
law. It is comprehensible that the probabilities of those events describe the properties of 
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random point sets in great detail since the bounded sets B in (1) can be chosen arbitrarily 
small. For calculations, it is very convenient to express a point configuration 77 as a sum of 
Dirac measures, 

00 

1=1 

Since 77 is assumed to be locally finite, one can always find a suitable sequence a^, i = 1, 2, . . . 
of points in IR*^, e.g., by collecting the points in centred balls of increasing radius and giving 
them consecutive labels. 

An interesting class of point processes are the stationary ones, where the probabilities of 
the events (1) are invariant under translations of the bounded sets B. The best explored class 
of point processes is the class of Poisson point processes and among them the stationary ones 
in particular: Let A G R'*' and let vol(-B) denote the volume of a Borel set B. The Poisson 
point process with intensity A, Px, assigns to our typical events the probabilities 

Pa ({?? e ^•(R'^) has k points in s}) = exp (-A vol(B)) . (2) 

The expected number of points in a unit cube thus equals A. It is clear that Px is stationary 
because the probabilities only depend on the translation invariant volume. 

In the general case of Poisson point processes, the volume in (2) is replaced by an arbitrary 
locally finite measure p on IR*^. This results in the Poisson point process with intensity measure 
p, denoted by Pp. The process is stationary as long as the intensity measure is translation 
invariant, which in the case of R"^ just means that Pp = Px for some given A > 0. For special 
intensities, it might happen, with positive probability, that point configurations have more 
than one point in one place. Such situations are excluded if the intensity measure p has no 
pure point part, i.e., if p({a}) = for all a G IR*^. Such point processes are called simple. 
While Px describes some ideal gas, one might interpret Pp as gas of non-interacting molecules 
in a certain physical potential. In contrast to general Gibbs measures (cf. [5]), the particles 
in the Poissonian case are always non-interacting. 

3. Clusters and Cluster Properties 

Similar to the above mentioned models for an ideal gas, randomness enters the approach 
of this paper by means of point point processes, as described in the previous section, where 
there are a lot of well-known constructions and simulations [21, 20]. But the information one 
gets out of the typical construction rules are of a more global nature, like distributions. To 
describe the local properties of the modelled objects and still not to lose the randomness of 
the point processes, we need a proper concept, which will be described in this section. 

As mentioned before, a cluster from our point of view is a finite subset of IR*^. Let 

Xd := {x C R'^ I card(X) < +00} 

be the space of clusters in IR'^ (here card(X) denotes the cardinality of a set X). Typical events 
in this space are constructed analogously to ^*(|R^). The method to consider (random) 
collections of clusters is inspired by the theory of random sets by Matheron [14] . 

To attach clusters X e Xd to point configurations 77 G ^'(IR'^), we use the concept of 
cluster properties. A cluster property ^ is a measurable (cf. [22] or [15]) STibset of the product 
space Xd X ^* (IR*^) . The elements {X, 77) G are the clusters and point configurations which 
are 'connected' in the context of the connection rule Although in general misleading, it 
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Figure 3. The clusters of type in a point configuration might be inter- 
preted as a collection of non intersecting balls 

might be helpful to imagine the (random) points 77 as a set of nuclei, and a connected cluster 
X as the orbiting electrons of one nucleus, where - of course - not all configurations are 
possible. 

If we take Voronoi and Dclone tessellations as typical objects we want to describe, it is 
easy to see that, for a given point configuration 77 G ^'(IR'^), we need two concepts for the 
clusters. In the case of the Voronoi tessellation the vertices of the cells are different from 
the generating point configuration, while the vertices of the cells in the Delone case coincide 
with the generating points. To give those concepts a name we will call X G Xd a duster 
(of Type Q)) for 77 if just [X^rj) G 3). If, additionally, X C r] we will call X a cluster in rj. 
In this notation, the vertices of the Voronoi cells are just clusters for the generating point 
configurations, while the vertices of the Delone cells are clusters in it. (To stay in the image 
of nuclei, clusters in a configuration might be imagined as neighboured or even interacting 
nuclei, where the cluster properties are the interaction rules.) 

It can also be convenient to extend the concept of cluster properties to a situation where the 
clusters and the point configuration do not lie in the same underlying space IR'^. Therefore we 
will also consider measurable subsets of Xd x ^'(IR") as cluster properties, where d^n. In 
this context, only clusters for a point configuration are well defined. We will see an example 
below. 

Before we get to more difficult examples and tilings, let us consider a simple cluster property 
for illustration. Let r G IR"*" be fixed. The pair [X^rj) is defined to belong to the cluster 
property C x ^* (IR^) if X consists of a single point ax and rj has no point which is 
closer to ax than r except possibly ax-, when this is also an element of rj. One cluster of type 
Sir for a given configuration is easily interpreted as a ball of radius r which does not intersect 
with the points of 77. But the collection of clusters for 77 has no reasonable interpretation. The 
collection is also not locally finite in the sense that only finitely many clusters intersect with 
arbitrary bounded sets of IR''. On the other hand, the collection of clusters in 77 is locally finite 
and can be interpreted as a collection of balls of radius r/2 around the points of 77 where all 
the intersecting balls are removed. Imagine 77 now as a random realization through some point 
process. Then, the collection of clusters in 77 'is' a random collection of non intersecting balls 
(see Figure 3). Note that in the general case, where clusters in or for a point configuration 
consist of more than one point, the clusters may intersect. 

If one wants to examine collections of clusters of certain type for or in random point 
configurations, it is interesting to know whether the underlying probability law, namely the 
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chosen point process, produces interesting collections of clusters. One indicator for this would 
be infinitely (but locally finitely) many clusters for almost all point configurations (with 
respect to the point process). In the case of the clusters in a random configuration, there is 
the following interesting result from [22]: 

0-oo-Law of Stochastic Geometry. Suppose that the cluster property is of the kind that 
translating a pair {X, rj) does not alter whether it belongs to the cluster property or not. Let 
also P be a stationary point process. Then, with probability one, we find either infinitely many 
clusters in rj or none. 

In the case of the Poisson point process and subject to a mild extra assumption (see 
[22]) one has an even stronger result: Assume that with positive probability there exists at 
least one cluster in a randomly realized point configuration. This is sufficient to almost 
always (with respect to P\) having infinitely many clusters in such a configuration. Again, 
it is convenient and possible to express collections, this time not of points but of clusters, by 
means of Dirac measures. Here, the collection (p^{r]) of clusters in 77 might be expressed by 

X is a cluster of type S> in rj 

= ^l^{X,n)5x, (3) 

XCt] 

where 1^ denotes the indicator function of the cluster property ^. 

Although we do not have the 0-oo-law in the case of clusters for a configuration, we will 
see an example where the collection of all the clusters for random ?7's is locally finite. If we 
combine this function tpr/ with a point process P (more precisely, we take the image of the 
point process under the transformation ip^), we have a probability measure if^{P) on cluster 
configurations (again, see [15] for details). If we see one cluster as a 'point' in Xd, the notation 
^'{Xd) for the (locally finite) cluster configurations is sensible. A probability measure on 
this spaces is called cluster process. 

We will now collect some information about a special family of clusters. 

4. Geometry 

The cluster configurations of interest for this article are certain discretizations of tilings. 

Therefore, we need to adopt some well-known concepts of tilings to this case. First we need 
to define discrete polytopes. Since we want to identify a convex polytope with its vertex set, 
we take the properties of vertices for the definition: A cluster X E Xa is a discrete polytope if, 
for all points a X, there exists some hyperplane H such that the intersection of H and the 
convex hull {X) of X consists only of the point a. (See Figure 4 for illustration.) It is quite 
obvious that there exists a one-to-one correspondence between convex and discrete polytopes. 
The convex polytope is retrieved from a discrete one by taking the convex hull. 

Similarly, discrete simplices are obtained. For the construction in the next section, we need 
a strong property of simplices that is well-known in the 2-dimensional case: Every triangle 
has a uniquely defined circumcircle. In higher dimensions, a full-dimensional simplex {X) 
has a uniquely defined circumball K{X) where the complete vertex set - the discrete simplex 
X - is contained in the border, the circumsphere S{X) (see Figure 4); we refer to [15] for a 
proof. 
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FlGURE 5. Circumball and circumsphere for a simplex (in this case a triangle) 

A collection /j, G ^'{Xa) of discrete polytopes is called a cluster tessellation if the collection 
of the convex hulls of the clusters form a locally finite face-to-face tiling. In our context 'face- 
to-face' does not necessarily mean that there are no holes in the tiling. It just states, that 
if two tiles intersect, they intersect in whole faces. The collection /i is called simplicial if all 
polytopes in it are simplices. It is called complete if the collection of convex hulls covers the 
whole of IR''. If a cluster process is concentrated on the set of cluster tessellations, it is called 
random cluster tessellation. 

5. Examples 

In this section, we give three examples for random cluster tessellations, constructed via 
cluster properties. The first one consists of clusters in a point configuration and the second 
one of clusters for a point configuration. In both cases, the underlying space for clusters 
and point configurations is the same, the cluster properties are subsets of x ^'(IR"'). In 
the third example, the point configurations lie in some higher dimensional space, where the 
underlying space of the clusters can be interpreted as an embedded subspace. 

First example: a special Delone tiling. The cluster property that generates our tiling is defined 
as follows: Let i? G be fixed. For a discrete simplex, let K{X) =: K{X) \ X denote the 
circumball of X, where the discrete simplex is removed. 

A tuple {X,ri) & Xd x ^'(IR"') belongs to the cluster property S>r if 
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(i) X is a c?-dimensional simplex, 

(ii) K{X) does not intersect with rj, and 

(iii) K(X) has a radius < R. 

Figures 6 and 7 illustrate clusters of type in a given point configuration ij. The first two 



'e G r} 



Figure 6. Illustration of clusters of type in a given r] 




e G 77 



K{X) 




K{X) 



(a) X is not a cluster in r), because the rar (b) X is a not cluster in r;, because there 
dius of the circurnball is to big is another point of the configuration in the 

circumball 



Figure 7. Illustration of clusters that fail to be of type &r in a given rj 

assumptions are based on ideas of Delone [4] and ensure that, for a given r] G ^•(IR'^), the 
cluster configuration 

X is a cluster of type 5?^ in r/ 

is face-to-face and simplicial. The assumption (iii) makes the cluster configuration locally 
finite and thus a tessellation. On the other hand it produces holes in the tessellation (more 
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precisely: in the union of the convex hulls of the clusters) when the points in t] are not dense 
enough. (Figure 8 gives a typical section out of such an incomplete tessellation.) Thus, we 




Figure 8. A typical section of a tessellation of the kind of Eq. (4) 
get the following results: 

Proposition 1. Any simple point process P generates a random tessellation in the form of 
V&r{P) 0,^ explained in Section 3. 

If we again consider the Poisson point process P\, we can apply the 0-oo-law of stochastic 
geometry to obtain: 

Proposition 2. P\- almost surely, 

(i) there are infinitely many clusters of type S>r in a realization rj, and 

(ii) there are holes in the tessellation. 

Here, (ii) holds due to the fact, that a realization of a Poisson point process almost always 
has arbitrarily big gaps somewhere between the points. This, see [15], could possibly lead to 
models for random holes in certain condensed matter. 

Second example: discrete Voronoi tilings. This example describes a construction of the well- 
known Poisson Voronoi tiling. In [15], a generalization, so-called random Laguerre tessel- 
lations, based on marked point processes and the theory provided by Schlottmann [19], is 
constructed. However, the idea of random tessellations constructed as clusters for random 
point configurations is better illustrated in this easier case, so we stick to it. 

To understand the generating cluster property, recall the definition of Voronoi cells for a 
given point configuration 77 G (K'^) . If a G rj, the Voronoi cell of a in r/ is given by 

Vj^{a) := {v eR'^l No other point of rj is closer to v than a,} . 

We call a the center of the Voronoi cell 14? (a). The definition is illustrated in Figure 9. If the 
convex hull of rj is the whole space IR*^, the Voronoi cell is a convex polytope. 
We can now define our cluster property for this example: {X, rj) G if 

(i) the convex hull of rj is IR*^ and 

(ii) X is the vertex set of a Voronoi cell in ij. 
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Figure 9. The Voronoi cell oi a in r) = 6a + Sa^ + + ■■ ■ 

In this case, it is easy to see that, for every rj G (IR*^) , there are no clusters of type in 
r). But as long as (i) holds, the collection 

X is a cluster of type for jj 

is a complete cluster tessellation (a proof can be found in [21]). See Figure 10 for an illus- 
tration. The Poisson point process Px produces point configurations of the kind (i) with 




Figure 10. X = {vi, . . . ,vq} is a cluster of type for the configuration rj 
since it is the vertex set of the Voronoi cell V^r?(a) 



probability one. Thus, we have: 

Proposition 3. ip^y(P\) is a complete random tessellation. 
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Third example: a random cut and project tiling. This example is based on the so caUed cut 
and project scheme, a method to obtain tiUngs from a (higher dimensional) lattice. The 
vertices or the centers, in the sense of Voronoi cells, of the tiles are projections of subsets of 
the lattice. A detailed description of the underlying theory can be found in [16]. A large group 
of deterministic tilings can be constructed this way, for instance the Penrose Tiling (cf. [3]) 
or the Ammann-Beenker Tiling (cf. [2]). We present a way to construct random tilings which 
are 'close' to the known deterministic ones, where 'close' will have two meanings: in the first 
one the random points still lie on the lattice but not every point of the lattice will appear. 
The second interpretation will produce one point for every lattice point, but the points might 
be randomly shifted within some given radius. 

Again, we stick to a simple example, a 1-dimensional tiling. The mechanisms for random- 
ness can easily be adapted to any tiling that can be obtained via the cut and project scheme. 
The deterministic case of our example is described in [1]. 

Consider the 2-dimensional lattice 

A:= ^(u + v , u — V V2j u, V integersj C . 

Let vr be the projection onto the first coordinate and tt* the one onto the second. Consider 
the set 



W 



eeU^ 



7r*(e) G 



1 1 

71' VtJ 

which we will call strip. The projection tt{A D W) forms the vertex set of the so-called silver- 
mean chain, which is a deterministic aperiodic tiling (again cf.[l]). Figure 11 illustrates this 
concept. 



w 



fe£ A 



7r(e) 



Figure 11. The projections of the points in yl PI form the vertex set of an 
aperiodic tiling 



The appropriate cluster property C Xi x^* (R^) can be defined as follows: {X, rj) G Sig 
if and only if 

(i) cardX = 2, 

(ii) X = {a = 7r(ei), h = 7r(e2)}, with ei ^ e2 ^ rj ^~\W and 
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(iii) the intersection of the open intervall (a, b) and 7r(r/ n W) is empty. 
In this 2-dimensional case, it is easy to see that, as long as Tr{r) n W) is a discrete point set, 

X is a cluster of type &s for V 

is a 1-dimensionaI tessellation where neighboured points form a cluster, respectively the vertex 
set of a tile. 

To embed the deterministic version of the tiling into point process theory, just define Pa to 
be the point process in which produces the lattice A with probability one. Then, ip^g{Py[) 
is a process that almost surely produces the silver means tiling. This tiling consists of two 

prototiles of length 1, and 1 + \/2, respectively. 

For the first random version, consider the discrete measure 

p:=Yc-Se, 

eeA 

where c is some positive constant. As mentioned above, Pp, the Poisson point process with 
intensity measure p, might produce point configurations with more than one atom in a single 
point, in this case in the points of A. All the point sets produced by the random mechanism 
Pp have the form 

ry = ^fce(r?) - de, 

where ke{r)) is some natural number or 0. The support of such an 77 is defined as 

supp(77) := ^ 5e. 

eeA, ke{ri)jtO 

The support is a subset of A, especially has only one atom at every point. Let P* be the 
image of Pp under the mapping supp. Thus: 

Proposition 4. P* is again a simple point process, where the realizations are random subsets 
of the lattice. The probability for a certain point e of A to be in the random set is 1 — e~'^. 

Figure 12 shows a typical randomly realized point set and what happens by taking the cut 
and project clusters. Since the holes in such a random subset cannot be controlled, the tiles 
corresponding to the clusters for a random 77 might have any length of the form n + m\/2, 
n, m non- negative integers. 

The next point process randomly shifts the points of the lattice. Another way to randomly 
shift the points is presented in [11]. There might be some way to transform these two ap- 
proaches into one another. For the construction of our point process, let e > and ^^^(e) the 
ball of radius e and centre e, e £ A. e should be choosen small enough so that the balls do 
not intersect. Consider the mapping be ■ ^*(IR^) ^'(IR^), 

I card(r,nB,(e)) E/6^nB.(e) / > V ^ B,{e) / , 
\e, if7?n5e(e) =0, 

which gives the barycentres of all the points in 77 n B^{e). The configuration 

6(r?) := J^beW 

eeA 
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Figure 12. P* produces random subsets of the lattice A, the tiles might 
become larger 

of all these barycentres, by construction, has exactly one point in every e-bah around the 
lattice points. Thus we have the fohowing result: 




Figure 13. For typical realisations r] of P5 the projections of rj CiW form 
slightly changed tiles or even completely new ones 



Proposition 5. // you take some arbitrary simple point process P, e.g. P\, the image P5 
of P under the mapping b randomly produces point configurations with exactly one point in 
every e-hall centred in the lattice points. 

Figure 13 illustrates the typical situation and the resulting projections. The corresponding 
tiles to the clusters for a given r] typically are close to the tiles of the original silver means 
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tiling, differing in length up to 2-£. But since certain e-balls of points in ylPlM^ intersect with 
the strip W, there might be 'completely new' tiles. Nevertheless, the density of the vertices 
stays the same, since the probabilities to shift a point into and to the outside of the strip are 
the same. 

In both of the cases of the third example of this article, it is easy to see the following: 

Proposition 6. The cluster processes ip^g{P*), respectively tp^g{Pb), are random 1-dimen- 
sional tessellations. 

6. Conclusions 

Point processes, especially the Poisson point processes, in combination with cluster prop- 
erties give access to modelling discrete random structures. The information of the cluster 
properties in this case carry the information about possible connections respectively interac- 
tions of the particles. The third type of example shows a way to slightly randomise aperiodic 
tilings. In this context, future applications to random tilings in the sense of Gummelt [7, 8] 
are of interest. If point processes could be constructed that almost surely produce more spe- 
cial configurations, e.g. Delone or FLC (cf. for instance [13]) sets, the presented methods 
might get closer to applications like glasses or foams. 
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